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Weak gravitational leasing is normally assumed to have only two principle effects: a magnification 
of a source and a distortion of the sources shape in the form of a shear. However, further distortions 
are actually present owing to changes in the gravitational field across the scale of the ray bundle 
of light propagating to us, resulting in the familiar arcs in lensed images. This is normally called 
the flexion, and is approximated by Taylor expanding the shear and magnification across the image 
plane. However, the physical origin of this effect arises from higher-order corrections in the geodesic 
deviation equation governing the gravitational force between neighbouring geodesics - so involves 
derivatives of the Riemann tensor. We show that integrating the second-order geodesic deviation 
equation results in a ‘Hessian map’ for gravitational lensing, which is a higher-order addition to 
the Jacobi map. We derive the general form of the Hessian map in an arbitrary spacetime paying 
particular attention to the separate effects of local Ricci versus non-local Weyl curvature. We then 
specialise to the case of a perturbed FLRW model, and give the general form of the Hessian for 
the first time. This has a host of new contributions which could in principle be used as tests for 
modified gravity. 


I. INTRODUCTION AND OVERVIEW 

Weak gravitational lensing is becoming an important cosmological probe. The usual weak gravitational lensing 
theory depicts that a lens induces a convergence (a spin 0 mode) and a shear (a spin 2 mode) to a source lying behind 
it. There are two ways to describe this for a given mass distribution: one is to calculate all the null geodesics converging 
at an observer, and examine the output. Another is to calculate the propagation of a geodesic deviation vector, 
and examine the invariant moments of the resulting image distortion. This results in the Jacobi map between 
a source and an image [DE]. The ‘weak lensing’ calculation route typically takes the second option, as it produces 
accurate results easily and intuitively. The computation of the convergence and shear is achieved from the geodesic 
deviation equation which is linear in the deviation vector:^ 

r + = (1) 

Yet progressively stronger lensing events produce arcs and other more complicated distortions which cannot be cap¬ 
tured by a simple convergence plus shear distortion. How can these be described within the weak lensing formalism? 

The geodesic deviation equation is of course linear by construction: the rhs of 0 should read 0(^,^)^, for the 
terms that are ignored. More complicated lensing events can therefore be described by examining this equation to 
higher order. Up to second-order in ^ we have the Bazanski equation [3H6]: 

1“++ sAVfeC'r = o(e, o' ■ ( 2 ) 

Naturally, derivatives of Riemann induce higher-order changes in the deviation vector. In this paper we extend the 
general weak lensing formalism to include all second-order corrections in the deviation vector. We aim to give the 
general solution to the Bazanski equation in an arbitrary spacetime, in terms of a ‘Hessian map’ TLabc^ which is the 
higher-order equivalent of the Jacobi map Jab'- 

Ca = JabC^ + TLabcC^C^ (3) 

where is an angle between neighbouring rays at the observer, and is the same at the source. The Jacobi map 
is a rank 2 tensor in the image plane and, as such, has two spin 0 (a convergence from the trace and a rotation from 
the anti-symmetric part) and spin 2 degree of freedom (the shear from the trace-free symmetric part). The Hessian 
map is a rank 3 tensor in the image plane, has a symmetry on two indices, so can be irreducibly decomposed into 


^ a, 6, c, • • • denote spacetime indices. A, B, C, ■ ■ ■ are tetrad indices in the screen space, fc or g as an index denotes projection of that 
index in the direction of k or is the Riemann tensor. A dot is a derivative along the null curve — full details below. 
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two vectors and a symmetric trace-free rank 3 tensor. In general, these correspond to two spin 1 modes, and a spin 3 
mode (each having two polarisations). 

Secondary lensing effects have been considered in the past, and generically go under the name ‘flexion’, describing 
classic arc shaped images, though there are other effects from secondary lensing [THIS]. In these works, the flexion 
has been derived by Taylor expanding the amplification matrix (which is usually the linearisation of the Jacobi map 
about an FLRW background) across the image plane. In the notation above, this implies 

T-Labc^ (cJ\A\B) ■ (4) 


This gives an approximation to the full Hessian, involving screen derivatives of the convergence and shear. These are 
in fact the leading contributions in most situations, because they capture the changing gravitational field across the 
image place, integrated along the line of sight. 

The leading contribution to T-Labc in n perturbed Minkowski background can be derived easily. The leading term 
is the one with the largest number of screen-space derivatives, and comes from the third term in § 


6 AC 










(5) 


where all terms without 3 screen space derivatives have been ignored. Integrating and projecting everything into the 
screen space then gives 

Uabc - x" V^VbVc r dx' (6) 

Jo X 

which is the 3rd derivative of the usual lensing potential. This is the same as can be derived from Q. However, as 
we shall see there are many more contributions to the Hessian than given by this approximation. 

In this paper we make the link between flexion and geodesic deviation for the first time. In doing so we derive 
the general form of the Hessian, valid in any spacetime, which contains many more subtle contributions beyond the 
approximation Q. We shall then specialise to the case of perturbed FLRW model. 


II. DESCRIPTION OF A NULL CURVE AND THE SCREEN SPACE 


Consider a light ray with tangent vector and affine parame¬ 
ter A on the past light cone which is a constant phase hypersurface, 
S = const : 

ka=WaS. (7) 

The tangent vector is null and geodesic: 

fcafc“=0, fc'’Vfc/c“=0, (8) 

and may be decomposed relative to an observer with 4-velocity u°‘ 
into parallel and orthogonal components: 

= i-ubk'^) (u“ - n“) = F; (w“ - n“), 

nan“ = I, = 0, if = —Ubk^. (9) 

Here is the unit direction vector of observation, and is along 
an incoming light ray on the past light cone of the observer. E is 
the photon energy measured by u°‘. Note that our choice of is 
opposite to the direction of photon propagation. 

The screen space is orthogonal to the light ray and to the ob¬ 
server 4-velocity, and the tensor 

^ab — gab “f Ua’^b BaBb , (10) 



Figure 1: Spacetime diagram showing our defini¬ 
tion of 4-velocity and null vector, x S'lid g are used 
in Sec. lYl 


projects into screen space. It satisfies the following relations 

N\ = 2, NacN^b = Nab, Nabk’^ = NabU^ = NabU^ = 0 . 


( 11 ) 
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Given the 4-velocity u“ we can invariantly decompose tensors into scalars, vectors and projected, symmetric and 
trace-free tensors [SUES], which are all spatial. Then, For any spatial tensor T“ '...b, we can isolate the parts lying 
in the screen space and parallel to n“ in a 1-I-1-I-2 decomposition [551151] : 

( 12 ) 

7), = na---n^T^-...b. (13) 

For a PSTF tensor, we use |j to represent all the indices which are projected along n“. Traces in the screen space can 
also be removed. For example, a rank 2 PSTF tensor has the full decomposition 

Xab = X/^ab) = ^11 [naTlb - ^Nab^J + 2X|| ■ (14) 

The invariant decomposition of the covariant derivative of the photon ray vector is given by, 

Xbka = ^dNab + (Tab , (15) 

where 

e = N-'^Vakb, (Tab = N^a^Nb)'^Vakd - \eNab ■ (16) 

Thus 9 describes the rate of expansion of the area of a bundle of light rays and aab describes its rate of shear (the 
trace-free part of the derivative projected into the screen space). Note that there is no null vorticity since ka = VqiS. 
The covariant derivative of the 4-velocity is invariantly decomposed as: 

XbUa — XaUb '^^{Sab 4” UaV’b') Xjah T Hab i (1^"^) 

where 0 is the volume expansion rate of the worldlines, Aa is the 4-acceleration, Y,ab is the shear tensor and flafa 
is the vorticity tensor. In terms of these variables, the 1-I-1-I-2 decomposition of the covariant derivatives of E and n“ 
are: 


E aE — —A||Uq-|-—© ria + S||^-|-f2||^ 

XbUa = -A|| {UaUb - UbUa) + ^©(Uanfe - HaUb) + {Ua - na)(S||^ -|- flyj,) - UbAa 

1 \ 9 \ 

+ ( o© ~ nV; I Xab — E ^(Tab + X,ab + ^ab 


^3 2E 

The Sachs propagation equations for the null shear and null expansion are found from the Ricci identities 

_ ab _ p h°‘lA 

2 ^ -^ab^ ^ i 


(18) 


(19) 


DcTab 

dA 


— ~9(Tab Cacbd.k'^k^ , 


( 20 ) 

( 21 ) 


where D/dA = k'^Xa- The photon geodesic equation ^ reduces to the equations for the photon energy and observa¬ 
tional direction 


^ = -F 
dA 


^ - Aan°'+ T,abn°-n^ 


Dn“ 

“dA" 


= E 


{Abu^ - ^bcTi^n^) - -k (E“b - fl“t,) n'’ -f ( -0 - Abu’^ + EbcU^E 


1 


( 22 ) 

(23) 


We shall project vectors and tensors onto the screen space using the tetrad basis A,B = 1,2 where 

e^’’ub = e^i^nb = 0, NabeA°'eB'^ = gabeA^'es^ = Sab, NabCA^ = 0. (24) 


The derivatives of the tetrad (Ricci rotation coefficients) are not required here. 
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III. THE NULL GEODESIC DEVIATION EQUATION TO SECOND-ORDER 


Consider a deviation vector lying in the screen space orthogonal to k°‘. This links two neighbouring geodesics, 
one at x{X) and the other at x'(A). The deviation vector at x{X) is defined as the covariant derivative of Synge’s world 
function, which is half the squared proper length of the unique geodesic connecting x and x' - see . To second-order 
this obeys a generalisation of the GDE, known as the Bazanski equation (' = D/dX): 


C + k’^Ck‘^R\cd + + 2i^^^k‘^R‘^bcd = 


(25) 


An additional term ^^°^'^^'^R‘^bcd would make this accurate in all powers of but we shall not include this term here. 
Our initial aim is to give the solution to this equation. 

We shall project this onto the screen space using the tetrad basis The tetrad components of in the screen 
space are = e^^b- This is called the Sachs basis. The full connecting vector is 

ia = + ^kka ■ 


(26) 


(27) 


This has a part parallel to k°‘ at second-order which cannot be set to zero as it can at first. This part obeys 

ik = -l^^^^RAcBdk^^k^ + 2i^eRAcBdk^k'^ . 

Thus, the deviation vector is forced out of the screen space as it is transported along k°'. However, since is second- 
order, it can only influence the screen parts of via the linear term in the GDE - but by the symmetries of the 

Riemann tensor, this contribution is zero. We therefore do not consider any further. 

Then, the screen-projected part of (251 can be written as 


^A = R-ab^^ + 'PABci^i'^ + Qabc^^^ 


tBcC 


where 


TZab = —RABck'^ = R-ba ■, 


1 • 


1 


RaBC = --aRbAC — -A^A°‘^{B^C)^ ciRadbek^^k^) 


1 3 

QaBC = -j^RABC — ^RbAC ■ 


Here we have defined the Riemann tensor with one index projected onto as 

R\d = R^\dkb 

= C^\^kb + 5\R\kb - - \R5\ka] , 

where Rab and R are the Weyl tensor, the Ricci tensor, and Ricci scalar. Then, 


(28) 

(29) 

(30) 

(31) 

(32) 

(33) 


Rab = —CAcBdk‘^k'^—-6ABRcdk'^k‘^-. 

RaBC = CAdBck'^ + ^A[BR^C]kb ■ 

Rab is often called the optical tidal matrix. 

Now, the first-order part of is Lie dragged along [S], so obeys 


ia = ^^^bka = ab^^ 


^A = + ctabC 


In matrix form, 


^ = 3$, where Sab = ^OSab + cfab = 


^6 + ai (72 
(72 ^0 - (7i 


(34) 

(35) 

(36) 

(37) 


where CTi = an = —( 722 , (72 = a^ = a 2 i- S is the (linear) optical deformation matrix. Gonsequently, we have the 
alternative form for the 2nd order GDE: 


^A = Rab^ +Tabc^ C 
1 1 


BfC 
f). 


where 


TabC — 2^iBC)A — 2^A°'^{B^C)"^ciRadbek k^) + ^0 R(^bC}A — ^(RADiB + ^R(B\AD\)f^ C) 
We shall call this the (optical) distortion tensor. 


(38) 

(39) 
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A. The general solution of the second-order GDE 


observer, so that we have initial conditions 


We now turn to the solution of (|38|). We are interested in the case where the geodesic congruence converges at the 

(40) 


Ca(Ao) = 0, ^a(Ao)7^0. 


Since (38) arises pertubatively in powers of we can solve it perturbatively by writing 

e(A) = 6(A) + 6(A), 

where 6 = ^>(6)^ in the usual way. The linear part is the solution of 

Ca = 

with initial conditions at the observer 

Ca(Ao) = 0, Ca(Ao) = Ca ■ 

The solution may be written in terms of the (linear) Jacobi map, J^ab 

Ca(A)= Jab{\)C^ 

where Jab satisfies 

Jab = T^acJ^b^ Jab{^o) — 0 , JABi^o) = —Sab ■ 


(41) 

(42) 

(43) 

(44) 

(45) 


The Jacobi map takes a deviation vector at the observer and maps it to the deviation vector at the source. In matrix 
form we have 


j = nj, 


(46) 


which we shall find convenient later. Now, the linear also obeys (37), because it is Lie dragged along the geodesic 

(47) 


congruence - this relates the derivatives of ^ to those of Writing this in terms of J we have 

j = SJ. 

This implies the optical deformation matrix obeys 

s + ss = n. 


(48) 


Thus the Jacobi map is found by integrating (46), and is related to the ray bundles expansion and shear through ( |48| ). 

We shall now find the quadratic part of the non-linear GDE, assuming we know the solution to the linear part. 
Inserting the linear solution in terms of the Jacobi map into (28), ^2 satisfies 


^A — T^abJ = Fa 


FabcJ^J^ + QabcJ^J'e 


B AC 


_ n- rrB rrC /-D/-E 

— lABCJ eQ C 


Da-E 


rc 


or in matrix form 


^-n^ = F 


(49) 


(50) 


We can solve this by a modified variation of parameters method. For this we need an independent solution to the 
homogeneous equation ^ — 7?.^ = 0 in addition to J. Define a reciprical Jacobi map K.: 


k = nK with K{\o) = I, k{\o) = 0 . 


(51) 


This matrix can be found in terms of J and S. First, we note that since 71. is symmetric, 7?. = 7?.^, K, and J must 
be related by 


J K-J^K = -I 


( 52 ) 
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which can be checked by differentiating, and by verifying the initial conditions. This implies that K, obeys 


K = SK + {J^)-^. 


(53) 


We shall assume this solution is known, in addition to J. Now we suppose there is a solution to pOl of the form 


i = Kcx + J(3, 


(54) 


where we shall assume 


KdL + J^ = Q. 


Then (50) becomes 


fCa + j^ = F. 


In matrix form we are then solving 

The inverse of a block matrix is given by 


'A B 

-1 

[A- BD-^C)-^ 

1 

1 

1 

0 ) 

1 

C D 


-D-^C{A- BD-^C)-^ 

{D - CA-^B)-'^ 


which implies 

=( (/c-jrj'-i/c)-i 


We only need the first column of this block matrix. In particular, 

{k. - = \SK + . 

Therefore we have 




The second-order solution is then 


I = ^ dA' [x;(A) - j^{k)F{x'). 


(55) 

(56) 

(57) 

(58) 

(59) 

(60) 

(61) 

(62) 


The full solution to the second-order GDE can be written in terms of a linear Jacobi map and a quadratic Hessian 


as 


= JabC^ + 'HabcC.^C^ 


(63) 


where the new Hessian part is 

nABc{x) = j\x'[ic/{x)-jj,’^{x){j-^)j,^{x')jc/{k)]j%{k)jG^{k)jH^{x')TiBc{k). 


(64) 


One can check by differentiating that 


a ABC = Ti-A^'HoBC + TaDeJb^^C^ > 


( 65 ) 


which gives the differential relation between the Hessian map and the optical distortion tensor. 
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B. Extended weak leasing formalism 


We shall now outline the general procedure for finding the Hessian map in terms of the usual weak lensing variables, 
the convergence and the weak lensing shear (as opposed to the shear of the geodesic ray bundle). 

The Jacobi map can be expanded in terms of a mean distance and an amplification matrix, 


JT” — — dj\ 


1 - K - 7l 72 

72 1 - K + 7i 


det JT” 
det A. 


A 


( 66 ) 


(We are now assuming that the Jacobi map is symmetric, but extending to the case with rotation is straightforward.) 
It is normally assumed that this decomposition is around an FLRW background, by identifying 


dA 


det JT” 
det A 


(67) 


as proportional to the background area distance (proportionality determined via ^^(Ao) = —1), but really this 
decomposition of JT” is completely general. Within this ‘weak lensing’ interpretation, k is the convergence, and ^ab 
is the trace-free shear, where 71 = —711 = -I- 722 , 72 = 712 = 721 - 
We can expand the reciprical Jacobi map K. in the same manner: 


K. = dAA = dA 


1 - K - 71 72 

72 1 - k -I- 7 i 


( 68 ) 


where we shall assume that (iA(A) is determined from JT in the case when k = 7 = 0. From (521, this implies 


dA 

dA 


1 

dA dAdA 


(69) 


with dyi(Ao) = 0. The background area distance and reciprical distance are determined from the background expansion 
rate 9 via 


cLa = 


d-A — 2^^-^ 


‘-1 


(70) 


With this, the reciprocal weak lensing variables kAab are determined from (53), which becomes the system of linear 
differential equations: 


K + ( 


\dAdA 2 
^ ‘ 


--A6>) k - (Ti% + 0-272 = 


1 — (1 — k)^ 


dAdA 


7i 


72 + 


\dAdA ^ / 

\ -- 0-1 ) 72 + a-2k - 0-271 = 0-2— 

\dAdA ^ J dAdA 


7iA^ 

dAdA 

72Ai 


(71) 

(72) 

(73) 


We have defined A6 — 6 — 6 for convenience - but there are no approximations made. 


The weal lensing convergence and shear are found in terms of the Sachs optical scalars from (47), which become 

(74) 

(75) 

(76) 


k - ^A 6 »/c - 0-171 J- 0-272 = 

71 - ^A 6 » 7 i - o-iK = - 0 - 1 , 


72 - ( J- 0-1 ) 72 -f O- 2 K - 0-271 = Cr2 ■ 


The Sachs optical scalars are found from (48). The trace of (48) is 

9+^9^ + 2iaf + aj) = -Rabk'^k'^ 


( 77 ) 


















while the trace-free part becomes 


^AB + Sctab = CABcd.k^^k'’' , 
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(78) 


which are just the Sachs optical equations as derived from the Ricci identities, but now derived from the GDE 
combined with the fact that the linear deviation vector is Lie dragged along the congruence. 

The Hessian of the transformation is an integrated projection of the optical distortion tensor, 

T-Labc{^) = 'Pa^\Tdbc]{^) or 'Hbc(A) ='P[T_bc](A) , (79) 

where we have used a matrix notation for the first index of the Hessian and distortion tensor. The projection operator 
V is 


T = 


^ dA' [X:(A) - J{X)J^\\')K{}^)\ 




dydAi^T dAWdA(y)A(A)A (y) - dAWdA(y)A(X)A-\A') A(y)Ayy) 


(80) 


We now have everything in place to construct the Hessian of the lensing map. The strategy is to construct /C and 
J' as follows: 

1. Solve for 0 and then cIa and cIa- 

2. Find the Sachs optical scalars 0 and (Tab- 

3. Find the components of the amplification matrix, k and jab- 

4. Find the components of the reciprical amplification matrix, k and jab- 


5. Calculate the components of the optical distortion tensor (96). 


6 . Construct the integrated projection operator P (80). 


7. Integrate the distortion tensor against the projection operator to calculate the Hessian from (79). 


C. Influence of the Ricci vs Weyl curvature 

We can expand the optical tidal and distortion tensors in terms of the general matter variables and the Weyl 
curvature tensor, which we shall split into its invariant electric and magnetic parts. The Ricci tensor and scalar are 
related to the energy momentum tensor via 

Rob = Tab — \Tgab + Ag^b (81) 

The stress energy tensor has the invariant decomposition given the observer velovity 

Tab — P'^a'^h T pk-ah 4“ ^ab ■ (82) 

Then the fluid part of the Riemann tensor is 

{RfiuidT^cd = f (p-f 3p- -f I (p-f A)/l“[c/l^d] 

- h^\cqd] - - 2u[“u[c7r'’ld] -f 2. (83) 

The Weyl tensor can be decomposed into an electric and magnetic part: 

Eab = Cacbd yX u‘^ Hab = 5 Sade C‘^‘'bc u'^ , (84) 

so that the Weyl tensor is expanded as 


= + • 
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The projections of these tensors we require are: 

Rabk°'k^ = KS^(/9 + p + 2g||+7r||) , (85) 

Rj^kb = -KE{qA+TTn^), (86) 

CAcBdk^k'^ = (^SabE, + 2Eab - 2sciAHsf) , (87) 

CAdBck'^ = E (2SA[cEiig^+ ( 88 ) 

so that the projected parts of Riemann which appear in the GDE are 

E ^TZab = ~2^ab (p E p + 2^11 + TTy + 2E^^) — 2Eab + ‘^^c(aHb)' j (^9) 

E~^Rabc = -SA[Bi<lc]+'^iic]+‘^^i\c]) (90) 

E ^Rabc = —'^yi[B(9c] +'^iiC] + 2-E||C]) + £bc^iia 

”*""^(3 ^”^"”*"^") + t^\\c]E 2E\\(j^) — £bcH\\j^ . (91) 

Here we are using compressed notation whereby 

E\\a = ~ ^ ['®ii^(^ii ~ ^ 11 ) “f (—A^ + Sy^ + D,^I‘)Ei,a] + n^EbA ■ (92) 

We also require the projected derivative of the optical tidal tensor: 


E ‘^eA°'eB^^c‘^'^o-iRbdcek'^k^) — 2 '^-BC'[Va(p + P + 2(j'y + TTy + 2Ey) + 2(I]||^ + r2||^)(p + p + 2(711 + TTy + 2i?||)] 

-2 + Esa(C - ^A{C^ ( 2 £^|IB) - ^^llB)) 

'^E '^ 11 S) 

+4(I]||^ + fly^) (Ebc + Hbc) + 2e^“eB'’ec" {vaEbc - £d(cV|,|77,)^) . (93) 

We have used an over-bar to denote the reverse parity in the screen space: 

Xa = £abX^ and Xab = Sc{aXb) ■ (94) 

Consequently, the optical distortion tensor may be written as sum of parts induced by Ricci curvature Tabc^ 
non-local Weyl curvature, Tabc^ where: 

E ^'Tabc = 2^Mb “ ^ (^C) + 2 Sy^^-I-( p- l-p-I- 2 gy-I- ttii)-I-- (gyc)-I- ttiic)) 

/ 2 \ 

( 3 ® ~ ~ ^ j (9||B + ^I|b) 

+^BC - ^ (9iiA + ttiia) + - ^11 + E|| - —^ (^11^ -I- TTy^) 

+ ^ 0 -A(B ('Zlic) + '^IIC)) “ 4 ‘^BC {q\\A + '^11 a) 

(l^SciB + CrC(B^ ( 91 IA) +’"llA)) (\&^B{C + CTBic) (^llA) + ^11 A)) (95) 
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Figure 2: The 6 independent degrees of freedom of the Hessian, as defined by separating the anti-symmetric and trace-free 
parts of the tensor. This is the action on a unit circle. 


E-^T, 


ABC 


= -5 


E\\c) ~ E c) + 


c) 


IIC) 




C), 


E,-E 



30 ^ 


1 

• 

^1^ 

30 \ 1 

- 9iT + E ^-0 - H|| -1- E|| 

~ Ae j 

^\\A 

+ 2^a{b 

-H\\c) + ^ 

y3^ “ ^11 + ^11 

-1 

o 

1 


A{B 


1 0 1 ' 

30<5(C(B) + S(C(B) - ^iC{B) - -^S^dB) - ^<^{C{B) £^|Ia) 

+ 2 ('^^(s^^iic-) “ Sctbc'-E^iia) “ + -^eoAcr^d^iiB) 




1 D 


+ 


3 ®'^(C(b) 


S(C(B) - ^{C{B) 


H, 


-2E (e^b + ^ii{b) {e:c)a + Hc)a) + Eej^°‘edeQ^''Vb(^-Eca + £d{cHJl'^ . 


A) 


(96) 


In the last term in (96) we have refrained from using a separate notation for the fully projected derivative to avoid 
possible confusion. The double symmetrisation means 2{C{B)A) = C{BA) + B{CA). 


D. The effect of the invariant parts of the Hessian 

In general, a rank 3 tensor can be invariantly decomposed into antisymmetric, and symmetric trace-free parts [32j . 
In the case of Habc, which is a rank-3 tensor in 2 dimensions satisfying Habc = 'Ea{bc)-i the invariant parts are: 
Babc, which is totally trace-free, and 2 vectors Ha and Ha, where 

Ha = -^S^^Habc 

'Ea = {'^l-[BC]D +H[b\D\C]) 

Then the Hessian may be written 

Habc = H^a^bc) + £a{bHc) + Habc ■ (99) 

This gives a total of 6 degrees of freedom - 2 in each component. The number of indices indicates the spin level of the 
invariant parts - there are two spin 1 parts, Ha and Ha and one spin 3 part Habc- These correspond to F flexion, 
twist and G-flexion in the language of |15j . Because there are traces in each term in EabC’ R-icci curvature can only 
induce spin 1 degrees of freedom, while (non-local) Weyl curvature can produce all 3, and is the only source of Habc- 
Given an the angle at the observer between neighbouring geodesics, the Jacobi and Hessian maps transform 
this at the source into 


(97) 

(98) 


^a = JabC^ + HabcC^C^ 

= (i(l - k)Ca + djABC^ + H(AdBC)C^C‘^ + sabHcC^C,^ + HabcC^C'^ - 


( 100 ) 

( 101 ) 
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From this expression it is clear the action of the Jacobi map: the trace induces a uniform change in area of a source, 
while the shear induces two area preserving elliptical distortions at 45° to each other. To illustrate how Habc, 'Ha 
and Ha change the image of a source, let us examine each in turn, and their action on a circular image. First, Ha- 


'H{a^bc)C^C'^ = {^A + ‘2'HbC,^ o) = gC^ 


Hx [2 + cos 26] + Hy sin 29 \ 
Hx sin 29 + Hy [2 — cos 29] ) 


( 102 ) 


where Qa = ICICa = |Cl(cos 0, sin 0) in Cartesian coordinates in the plane. This produces a shift in the position of the 
image, and a change in area. There is a compression along the direction of Ha, leading to a cusp for large distortions. 
Next, 


sabHcC^C^ = C^sabC^ HcC^ = 


HxSin29+ Hy[l — cos29] \ 
-Hx[l + cos29]-HySin29 ) ' 


(103) 


This produces a sagging like distortion, perpendicular to Ha, with the normal modes at 90°. Finally the purely 
trace-free part of the Hessian also has two independent degrees of freedom, Hxxx and Hyyy (all other components are 
proportional to one or the other of these since Habc is trace-free): 


HabcC^C^ = C 


'Hxxx cos 29 'Hyyy sfu 29 \ 
'Hxxx sin 29 'Hyyy cos 29 J 


(104) 


This produces a triangular distortion at low amplitude, with polarisations at 30°. These curves are all trochoids of 
higher complexity than the circle and ellipse which appear at linear order - see Fig. (A trochoid is the locus of 
points traced out by a point attached at some distance from a circle, which itself is rolling around a larger circle.) For 
large amplitudes, we can see cusps appearing representing the formation of caustics. 


IV. PERTURBATIONS ABOUT AN FLRW BACKGROUND 


We shall now linearise the Hessian around a flat FLRW model. We shall write our perturbations with respect to 
the Poisson gauge, where 


ds^ = [ — (1 -I- 2^)drj^ + (1 ~ 24')7ijda;Ma;^] , 


(105) 


where 7 ij = 5ij if x* are cartesian (z, j, k,... denote spatial indices, usually on the conformal Minkowski background 
when they are raised and lowered with jij). The observer 4-velocity u°‘ is perturbed as 


u 


0 


-(1 — 4>), u* = -V*u, Ui = a\/iV 
a a 


where v is the first-order scalar velocity potential which obeys 


(106) 


2a 


V = — 




-{'S' + H^) 


'i'' + H^ 
H' -H^' 


(107) 


where H = a! j a is the conformal Hubble rate, and we also use H — Qj?, = Hlain the background, and the background 
density is p = 2{H^ — H')/a^. We also have for the time and radial parts of the vector n“. 


SrP = Su^ = -dyV, Sn^ = 4', 

a ^ 

and similarly for k°‘ = E{u°‘ — n“), 

= {l + z)5E-{l + zf{<^ + d^v), 
Sk^ = -{l + z)6E-{l + zf{'i>-d^v) 


(108) 


(109) 

( 110 ) 


where the perturbed energy is given by (enforcing that the time component of the photon 4-vector is unperturbed at 
the observer) 


6E — 6z J- (1 -l- z)($o -|- d^Vo) 


( 111 ) 
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with the redshift perturbation given by 


5z 

1 + z 


(d^v 


Jo 


Note that = 0 implies Sk'^ + Sk^ 

The density contrast 6 is given by 


-(1 + giving (1101. 


( 112 ) 


a^p6 = 2V2^' - 6n{'S' + n<^). 


(113) 


For a perfect fluid in GR we have 'I' = <&. However, we will keep these potentials separate for generality. In that case 
we have an effective pressure and anisotropic pressure perturbation given by 


a^6p = 2^'" + 2'H(24'' + $') + 2(2'H'+'H^)$-^V2(^'-$), 
TTij = V(,Vj)(4 '-$). 

In the LCDM case in GR, the first leads to the Bardeen equation $" + SHd)' + a^A$ = 0. 
We also require the matter shear and electric Weyl tensor. These are 


a 2 

Sy = aV(iVj)U = aVjVjU --7ijV u, 

= ^V(,V,)(vk + $) = iv,V,(vI/ + <i>)-^7,,V2(d/ + ci>). 
We denote the radial co-moving coordinate as y, where 


an^Wi 


A 

dx 


A + 

dx dr] 


(114) 

(115) 


(116) 

(117) 

(118) 


in the background. Because Tabc is first-order we do not need the perturbation of this. The derivative d/dy is along 
the null geodesic, so that x is a conformal affine parameter from the observer to the source obeying dA/dx = —a^. 
Together with this we have the spatial Laplacian, 

= 0-^7*^ = a-^V^ , (119) 

with the notation such that is only used to denote the co-moving spatial Laplacian (i.e., the Laplacian on a 
Minkowski background). This is expanded using Nab into a co-moving screen Laplacian and radial parts, as 




- Nl+dl + -d^ 

\ 

= A + ” 1 — 

dx dx X 



( 120 ) 


We also require a comoving tetrad in the screen-space (i.e., with the scale-factor factored out), which is the natural 
tetrad for the observer. We define these as 


06 /= a{n, e^} = 6/= {e^, Eyi} where / = X)^, and A = ‘d,ip. 
With respect to a standard Gartesian basis i,j, k, we have 


and 


6 ^ = sin 'd cos pi + sin § sin p j cos i) k , 
6 ^ = cos iJ cos pi + cos a sin pj — sin -d k , 
Bip = — sin pi -\- cos p j , 


V = e^d^ + —d^ + 


X 


_ 1 _ d 

XSini7 


( 121 ) 


( 122 ) 

(123) 

(124) 


(125) 
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The Ricci rotation coefficients for this tetrad are defined by 
V/ej = (e/ • V)ej = jjBk =» 


1 


pA pX ^ 

^ ^ AA- ^ 

P =-P ^_coP 


A 

r^BC" = ‘^sabSq 


cot d 


(126) 


X 

(no sum implied in the expressions on the right). The tetrad alternating tensor is eAB = ^b^- using 

notation where V/ = e^Vf,, a derivative which we only use on scalars or tensors projected onto the tetrad basis. In 
this basis the 3-d Laplacian is 

V/V^ = ■ V)(ej • V) - ■ Vej) -V = dl+ (127) 


which implies that the 2-D Laplacian is, 


wl = VaV^ - -d^ = VaV^ - - — + a, 


X 


X Vdx 


which in turn implies the 3d Laplacian is 


(128) 


= VaV'^ 


dx 


Radial and angular derivatives commute when acting on scalars as 

Va9p^ = 9^VaH—Va = —VaH—Va + Va9^, 
X dx X 

w d d ^ 1 

Vaa — A VaH—Va, 
dx dx X 

^ d 2 / d 2 \ d ^ 

'^dx^ Vdx'^xJdx 

The integral version of the second is found from xVAd/dx/ = d/dx(xVA/) as 


xVa [ dx' = / dx'x'VA • 
Jo Jo 


Projecting derivatives on scalars we have: 


1 


^A°'^B '^a^b = VaVb H-(5bai9x — P^^Vc , 

X 

Q 

h- CV7 V7 VT _ V7 V7_V7_ i Jv V7 ^ T^D 


= V^V^VcH- ^{AB^C)dx~^ BA^C^ D — C{A^ B)^ D 

X 


2T 


D 


E{C^ B)A 


- VaP^cb 


2 2 

Vb - --^(B\A\C)dx - -^Sa(b'^c) 
X X 


(129) 

(130) 

(131) 

(132) 

(133) 

(134) 

(135) 


Relative to the matter frame, the matter content is dust. For full generality we shall keep the effective pressure 
components non-zero so that our results can be used in the case of modified gravity. To linear order the magnetic 
Weyl tensor is zero for scalar modes. So, the distortion tensor becomes, relative to the co-moving matter frame 


Tabc — + 2)^'5a(b 


4 

“^(1 + z)^A{B 


(yc)P + 


IIC) 


(1 -I- z)Vc){p H- TTy) - T^wc) + ( “2^ + (1 + ) T^wc) 


-^Sbc 


^11^ ~ ( ~2^ + (1 + ) TTyA 


_Weyl 


BbC — + ^)^A(B \ + z)Vc)Ei\ + E\\C) + 


C) 


+ 2 ( 1 +-) 


-E^a H- ( H- 2(1 -I- z)E[ J i?||^ 


Sbc - (1 + zfeA’^e^^B^ecfVbE,. 


(136) 


(137) 
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The projected parts of Eat and Sab we require are: 

1 


2a^E, = 2e^e^Eab = 


2a^E^^^ = 2eXe°Eab = 


2a^e^ e^ge^-^WbEca — 


92 - -V 2 


($ + 4') = 


-sV.v^ 


3 Vdx 


($ + 4'), 


I “j—+ V/i(9^ 

[dx 


($ + 4-), 


(138) 

(139) 


V(bVc)Va — ”1- A)C>x ~ (BC)^D 

X 


2E^ A(^B^ C)^ D + {BC) + E(C^‘^\B\A) ~ ^ (B^^\A\C) 


~'Z) t^E 


D 




Vd 


--^A{BC)dx - -^^BC'^A - ■^Sa(b'^C) 
X XX, 


($ + 4-), 


2a^E«A — 


dx 


2% —va + - —va - — v^a 


dx^ 


dx 


(4> + 4-) 


VSm = -2(^^+a,) V^(4/' + H4>) , 
( 4 / - 4>) = 


a TTii = 


- 3^' 


-3V.V^ 




(4--$), 


dx 


+ v^a„ 


(4'-$), 


(140) 

(141) 

(142) 

(143) 

(144) 


Note the notation here: on the left hand side we are using the usual tetrad ej, while on the right the indices are with 
respect to e/. Really we should introduce new notation for objects in the tetrad e/, but we don’t need to: if an object 
is a 1+3 covariant tensor used above we’re using e/; if there’s 4>’s then the indices represent the comoving tetrad e/. 
We also require 


—VaP = + 


dx^ 


+ ai( + 2 (a„ H— 1 -— 


xj dx Vx 


--3-Ha 


Va4'-3H^Va$. 


(145) 


We can now give the perturbed version of the distortion tensor in terms of the metric potentials. Because this is to 
be integrated down the past lightcone, we write it in terms of d/dy before other derivatives using the commutation 
relations. The Ricci contribution to the Hessian becomes [indices on the left are with respect to the full basis, on the 
right the conformal one]: 


{l + z)-^7^lc = -;^<5A(B^Vc)VDV^4/-H(2^+2a^ + 3H) Vc)$ + 


dx 


-i^5a{b'^c) [d'" + 'a(24'' + $') + (2H' + H^)^] 


^-a^+2A 

dx^ X dx 


— — ^yl(BVc)V£)V +5a(B 


^ + -4H + - + aj — + -4H + - a 


dx2 ■ V ' X ' '^Vdx^' 


X 


— 3HI a 

X 


Vc) 


Vc)4' 


+5bc 


L dx^ 


-di^+ 4H---ajdT+ 4H-- a 


dx 


VaK^-^)- 


( 146 ) 


We have presented this such that in the case of LCDM and GR, only the first line remains. The Weyl contribution 
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becomes 


(1 + 2) ^Tabc = 2 1 “ + 2'^a(bVc)VdV^ + + 2r'°^(3Vc)V_D 

^a(b'^c) 


1^ 
2dx^ 
1 d2 


4^ 4x 2"^^ dx 


4x 


-^— I 9n H- 


X"J 


2dx2 


7, 


7 


1 


4x 


- — - -d^ ] — - -n - — ] ^ 


dx 


4x 


X^J 




Sbc'^ A 

E(C^ \B\A) ~ (Bi^^ \A\C) 


t^D t^E I "p-D t^E T-7 

EA^ {BC) ■*- iRM'i ~ V ('Ri 




(147) 


As the final part of the solution we construct the integral projection operator. This is simple because Tabc is 
already first order, so we just require it in the background where it reduces to the usual leasing kernel with small 
modifications. The area distance in the background is given in terms of the comoving distance 


dA{x) = 


X 


1 


dz' 


1 + 2 (1 + 2) Jo (1 + 2^)'H(+) 

from which we derive the reciprocal distance 

r 

dAM = dyi(A) lim 
£-> 0 + 


0 = --ii + zf{i-nx) 


dX' 


dA{Xo-e) Jx^_,dA{X'f 


X 


lim 

1 + Z e->0+ 

The projection operator simplifies to in FLRW, 


rx 1 

dx'^ 

£ X . 


1 


1 + z 


V = I 


dX'dAiXT \dA{X)dA{X') - dA{X)dAiX'] 


1 + z, 


dx' 


./ x' ix-x') 
il + z{x'W 


Finally we can give the full expression for the Hessian in the case of perturbed FLRW: 


'Habc = [ dx'x'^ix - x') 

^ z Jq 


C4Bc(x') + r^2c(x') 


(l + z(x'))® 

Note that the redshift terms cancel inside the integral with all integrated terms conformal comoving ones. 


(148) 

(149) 

(150) 

(151) 

(152) 

(153) 


A. Dominant contribution 

The dominant contribution to the distortion tensor arrises from terms with the highest number of screen-space 
derivatives, which are enhanced compared to radial derivatives on small scales. In this limit it is straightforward to 
write down the Hessian: 


1 rx 

= dxY"(x'-x)[V(BVc)VA($-bvl/)] (154) 

= ^(^^V(bVc)VaV’(x) (156) 
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The invariant parts for the dominant part now become 

Ha = 0 (158) 

Habc = ^(i^^(-4VbVc)V'- (159) 

Angled brackets in the last line denote the trace-free part. The spin 1 modes associated with rotation are not excited 
in this approximation, while the others are sourced directly by the density contrast. The spin 3 mode is sourced 
purely by the distortion of the electric Weyl curvature. 


B. Modified Gravity 


In principle, measuring the Hessian allows us to reconstruct and separately. Equality of these two potentials 
is a clear test of GR. Crucially, the pure trace-free mode of the Hessian Habc depends only on $ -|- tp, and is only 
induced by the Weyl distortion. However, the other two modes depend on $ -|- $ but also on and $ in different 
ways. The Weyl contribution always gives $ -I- d>, but the Ricci part sources the invariant parts of as: 


1(1 + 


= --VaVbV^($ + ^) 


11 d^ 
'TdP 
9 d^ 

8dx2 


while the rotational invariant part becomes 




4x 

21 


dx 


4x7 




+ ( — SH — 1— -dr) I — 


3(1 + 2) — + 


3^ 
'2dx2 ^ 
23 d2 
8 dx^ 


-7H 


4x 


-^ 7^1? 

X 4 


dx 


27 \ 


d 

dx 


~H + ^ ) d„ + 2H'+^H^ 

2 4x7 

17 \ 


V^$, (160) 


-¥^^4xJ 


dr, 


Wa'S 


+877 - - Idr, 

4x 4 


d 

dx 




:H^ 


(161) 


Consequently, careful comparison of the 3 independent modes can give new tests of modified gravity. 


V. CONCLUSIONS 

We have derived for the first time the general solution to the Bazanski equation, which is the extension to the 
geodesic deviation equation which underlies weak lensing flexion. Consequently we have given the general equations 
for flexion in a perturbed FLRW model. An additional mode is found which depends on radial derivatives of the 
potentials along the line of sight. It may be used to test for modifications of GR, as it depends on the potentials 
differently from the other two flexion modes. 
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